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Abstract 

We use covariant methods to analyse the nonlinear evolution of self-gravitating, non- 
relativistic media. The formalism is first applied to imperfect fluids, aiming at the kinematic 
effects of viscosity, before extended to inhomogeneous magnetised environments. The nonlin- 
ear electrodynamic formulae are derived and successively applied to electrically resistive and 
to highly conductive fluids. By nature, the covariant equations isolate the magnetic effects 
on the kinematics and the dynamics of the medium, combining mathematical transparency 
and physical clarity. Employing the Newtonian analogue of the relativistic 1+3 covariant 
treatment, also facilitates the direct comparison with the earlier relativistic studies and helps 
to identify the differences in an unambiguous way. The purpose of this work is to set the 
framework and take a first step towards the detailed analytical study of complex nonlinear 
systems, like non-relativistic astrophysical plasmas and collapsing protogalactic clouds. 

1 Introduction 

General relativity is believed to describe strong gravitational fields and also to determine the 
large-scale dynamics of our universe. Nevertheless, when the gravitational field is weak and on 
scales well inside the Hubble length, Newtonian gravity remains a very good approximation. 
The same is also true when dealing with low temperature (cold) plasmas, where the effects of 
special relativity are negligible. All these mean that Newtonian physics remains a very depend- 
able mathematical tool for a variety of astrophysical and cosmological studies. In particular, 
the theory can offer very useful insights regarding the behaviour of complex nonlinear systems, 
like a collapsing protogalactic cloud for example. Moreover, despite the fundamental differences 
between Newtonian and relativistic fluid dynamics, the two theories still share many close paral- 
lels. These analogies become more prominent and clear when using relative-motion descriptions, 
such as those the relativistic 1+3 covariant formalism and its Newtonian counterpart are based 
upon. Here we will use the latter. 

The covariant approach to fluid dynamics assumes the existence of a unique vector field that 
represents the average velocity of the matter at each point in space, or at each spacetime event 
in the case of a relativistic study. The formalism offers a Lagrangian description, where every 
kinematic and dynamic quantity is decomposed down to its irreducible parts; a splitting that 
combines mathematical compactness and clarity with physical transparency. The fluid kinemat- 
ics, in particular, are monitored through a scalar, a vector and a tensor field that respectively 
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describe the average volume evolution, the rotational behaviour and the shear deformation of any 
given fluid element. The evolution of these variables is determined by a set of three propagation 
equations, supplemented by an equal number of constraints. Once the full (nonlinear) expres- 
sions have been obtained, the covariant formulae can be applied to any physical environment by 
simply adjusting the symmetries. 

In the present article, we review the covariant approach to Newtonian hydrodynamics and 
provide the complete set of the nonlinear propagation and constraint equations that describe a 
bound, self-gravitating medium. We first consider the case of a barotropic fluid and examine the 
kinematic implications of inhomogeneity. This means looking at the gravitational collapse, the 
shear anisotropy and the rotational behaviour of the fluid. Our results show that overdensities 
tend to enhance the collapse, while underdensities act against contraction - or tend to accelerate 
the expansion. We also find that, under the barotropic-fluid assumption, vorticity cannot be 
generated. At each step, we compare our Newtonian expressions to their relativistic counterparts 
and establish the main analogies and differences between the two. Exploiting the advantages of 
the covariant expressions, we apply our nonlinear formulae to the case of an imperfect medium, 
in an attempt to investigate the role of viscosity. Among others, we find that a viscous fluid will 
generally act as a source of rotation. Also, by involving the internal properties of the fluid, we 
discuss how hydrodynamic flows can be represented as purely gravitational motions and outline 
the potential applications of this dynamical correspondence. 

With the full hydrodynamic equations in hand, we proceed to incorporate magnetic fields into 
our study. Introducing an electron-ion system and assuming overall charge neutrality, we derive 
the covariant magnetohydrodynamic (MHD) formulae for an electrically resistive fluid. These 
include, for the first time in the Newtonian limit, the covariant form of Maxwell's equations 
and allow for a direct comparison with their relativistic analogues. The evolution of the 5-field 
is monitored by looking at both the isotropic and the anisotropic components of the magnetic 
pressure. Confining to a barotropic medium, we consider the effects of the field on the fluid 
kinematics. The magnetic implications for gravitational collapse, for example, are encoded in 
Raychaudhuri's equation. The latter reveals how increases in the pressure of the .B-field assist 
the contraction by adding to the gravitational attraction of the matter. We also identify in 
covariant terms what is commonly referred to as "magnetic braking", and show how the effect 
results from the elasticity (i.e. the tension) of the magnetic forcelines. As with hydrodynamics, 
we take every opportunity to compare our Newtonian expressions to their general relativistic 
partners and identify all parallels and differences between the two sets. Thus, in contrast with 
general relativity, we find that the magnetic pressure has no effect on Newtonian vorticity. In 
agreement with the relativistic analysis, on the other hand, the magnetic tension is found to 
affect rotation and act as a source of it. Finally, by assuming a perfectly conductive medium, 
we apply our results to the ideal MHD case, establish the pattern of the magnetic evolution in 
such an environment and also discuss how the magnetohydrodynamic equations can be reduced 
to pure hydrodynamic ones. 

The main aim of this work is to introduce the key features of a formalism that will be 
subsequently used in nonlinear Newtonian hydrodynamic and magnetohydrodynamic studies. 
For this reason we have gone beyond the perfect-fluid approximation and incorporated viscosity 
effects into our equations. Similarly, the MHD formalism has been extended to allow for media 
of finite (nonzero) electrical resistivity. Our targets are nonlinear systems that are adequately 
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described by the Newtonian theory. These include non-relativistic astrophysical plasmas and 
protogalactic clouds (of subhorizon size) that have decoupled from the background expansion 
and started to collapse. 

2 Covariant hydrodynamics 

The covariant approach to fluid dynamics dates back to the 1950s and the work of Heckmann, 
Schiicking and Raychaudhuri [TJ. The formalism was originally applied within the Newtonian 
framework before extended to general relativistic hydrodynamics and magnetohydrodynamics 
(see [21 [3] for recent review articles and further references). In the present section we first 
review and later (see § I2.3ti2.5]) extend parts of where the reader is referred for more details. 
Relative to that article, there are also several notational differences, reflecting the presentation 
changes that have taken place since the early 1970s. For an alternative, 4-dimensional covariant 
approach to Newtonian hydrodynamics we refer the reader to [5]. 

2.1 Self-gravitating fluids 

We use fixed space coordinates {x a , a = 1, 2, 3} to define the metric tensor h a b of the Euclidean 
space, so that v 2 = h ao v a v b for any vector v a . The above given metric and its inverse h ab - with 
h ac h cb = 5 a b and 5 a a = 3, where 5 a b is the Kronecker symbol - are used to raise and lower the 
tensor indices. When one uses Cartesian coordinates, as we will be doing here, = 5 a b- Then, 
covariant and contravariant components coincide and partial derivatives are the 'correct' spatial 
derivatives (see [U [6] for further details) Q 

We adopt the fluid description, assuming the existence of a unique vector field representing 
the average velocity of the matter at each point. The 3-velocity field v a is tangent to the flow 
lines of the comoving (fundamental) observers. The time derivative of a tensorial quantity T is 
given by the convective derivative T = dtT + v a d a T, where d a = d/dx a . Thus, the convective 
derivative of the fluid velocity is 

v a = d t v a + v b d b v a , (1) 

with db v a describing the spatial variations of the velocity field (e.g. see [7]). Note that we adopt 
the Einstein summation convection, according to which repeated indices are summed. Like any 
second-rank tensor, the spatial derivative of v a decomposes as 

ObV a = g ®$ab + Vab + ^ab , (2) 

where = d a v a , o a b = d(b v a) an d ^>ab = <%,t> a ]ll The tensor d v a monitors the relative motion 
between two neighbouriong fluid flow-lineso In particular, determines the volume expansion, 
(Tab the shear deformation and uj a b the rotational behaviour of a given fluid element. Positive 

In a general frame h a b 7^ 8 a b and the covariant and contravariant tensor components do not always coincide. 
Then we need covariant, instead of partial, derivatives to compensate for the "curvature" of the system [4JE]. 

2 Round brackets in the indices denote symmetrisation, square indicate antisymmetrisation and angled ones 
define the symmetric and trace-free part of second-rank tensors. Therefore, d^Va) = dn,v a ) — (d c v c /3)S a b- 

3 The relative velocity vector (x a ), between two neighbouring flow lines, is related to their connecting vector 
(x a - connecting the same two particles at all times) via the transformation x a = x b dt,v a (see [1] for details) . 
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values for B correspond to an expanding fluid, while negative ones indicate contraction. The 
volume scalar can also be used to define a representative length scale (a) along the flow lines 
by means of a/ a = 0/3. In cosmological studies, the aforementioned length scale corresponds 
to the scale factor of the universe. The antisymmetry of the vorticity tensor implies that we 
can define a vorticity vector by means of ui a = e a b c ui bc /2, with e a b c representing the alternating 
tensor of the Euclidean spaceo By construction = e a i, c uj c , ensuring that oo a b^ b = 0. The 
vorticity vector also determines the rotation axis of the matter, namely the only direction that 
remains unaffected by the rotational motion [3] . Finally, the shear and vorticity magnitudes are 
defined by a 2 = a a ba ab /2 and to 2 = uj a b^ ab '/2 = ui a u a respectively [2]. 

Assuming that $ is the Newtonian gravitational potential, we use the velocity of the fluid 
to define the vector 

A a = v a + d a <$> , (3) 

which describes the combined action of gravitational and inertial forces. The vector A a corre- 
sponds precisely to the relativistic 4-acceleration and vanishes when the matter moves under 
inertial and gravitational forces alone [H [6]. The gravitational field is determined through a 
Poisson-like equation of the form 

d 2 $ = ^ K p-A, (4) 

where d 2 = d a d a is the Laplacian operator, k = 8ttG represents the gravitational constant, p is 
the density of the matter and we have allowed for a nonzero cosmological constant A (in units 
of inverse-time squared). 

2.2 Nonlinear hydrodynamics 

Using the convective derivative operator, decomposition ([2]) and definition ([3|), the nonlinear 
continuity equation and the Navier-Stokes formula associated with a self-gravitating fluid assume 
the covariant forms 

p = -6/9 and pA a = -d a p - d b ir ab , (5) 

respectively [I]. Note that p is the isotropic and tt^ is the anisotropic pressure of the medium 
(with 7r a b = TT( a b))- To close the system one requires the equations of state for the matter. These 
usually take the simple barotropic form adopted in § 12.31 or the phenomenological shape of 
Eq. (|22p in § 12.41 though in general they depend on additional thermodynamic variables. We 
also need a set of nonlinear formulae to describe the fluid kinematics. These comprise two sets 
of three propagation and constraint equations, which (like their relativistic counterparts) are 
obtained by applying the Newtonian analogues of the Ricci identities to the velocity vector of 
the fluid, namely by means of 

d[ t d b ]V a = and d [c d b] v a = . (6) 

4 By construction the volume element (the Levi-Civita tensor) has e a bc = £[ a bc], with £123 = 1. Also, e a bc£ dqp = 
Z\5 [a d & b q 5c\ v , which ensures that e abc e dqc = 2\S [a d S b] q , e abc e dbc = 2S a d and e abc e abc = 6. 
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The former of these integrability conditions leads to the propagation formulae. To be precise, 
the gradient of ([3]) together with definition ([1]) and decomposition ([2]), gives 

1 2 

{d b v a )' = -- Q 2 5 ab - - @(a ab + u) ab ) - d b d a ® + d b A a - a ca a c b + oj ca u° b - 2a c[a uj c b] . (7) 

This expression contains collective information about the kinematical behaviour of the fluid. We 
decode this information by isolating the trace, the symmetric trace-free and the antisymmetric 
components of (J7J). 

We begin with the trace of Eq. ([7]), which by means of leads to the Newtonian version of 
the familiar Raychaudhuri equation, 

<d = -\e 2 -l K p + d a A a -2(a 2 -u; 2 ) + A, (8) 

that determines the expansion (or contraction) rate of the fluid. Comparing the above to its 
relativistic counterpart (e.g. see Eq. (1.3.3) in [3]), we notice that only the density of the matter 
contributes to the gravitational mass and also note the absence of an A a A a -term in the right- 
hand side of (JSJ). 

In an analogous way, the symmetric and trace-free component of ([7]) provides the evolution 
formula of the shear 

2 

<jab = ~- &CF a b - E ab + d(^ a A b ) - CT c ( a Cr C fe ) + U c ( a UJ C b ) . (9) 

Here, E ab = d(hd a \$> represents the tidal part of the gravitational field and corresponds to the 
electric Weyl component of the relativistic treatment (compare the above to expression (1.3.4) 
in [3])H The vorticity term, on the other hand, carries the distorting effect of the centrifugal 
forces Also note that, in contrast with the relativistic analysis, there are no A^ a A b ^ and 7r ab 
terms in the right-hand side of (|9]). 

We close the set of the propagation formulae with the skew part of ([7]). The latter governs 
the rotational behaviour of the fluid element, either in terms of u ab 

2 

Wat = -g @u ab + d [b A a] - 2a c[a uj c b] , (10) 

or in terms of to a 

2 1 , 

u a = ~ g ©w a -- curL4 a + a ab LO , (11) 

since u> ab = e abc uj c by definition and cuvlv a = e abc d b v c for any vector v a . Therefore, when only 
gravitational/inertial forces are present and in the absence of shear, expressions (|10p . (jllh imply 
tj a oc a~ 2 and consequently constant angular momentum. Note that both of the above have the 
form of their relativistic counterparts (e.g. compare (|lip to Eq. (1.3.5) in [3]). 

Expressions ([8|)- (|lip monitor the nonlinear evolution of the irreducible kinematical quanti- 
ties of a Newtonian self-gravitating fluid in fully covariant terms. For a complete kinematical 

5 The tidal field can be associated with a component of the gravitational potential that does not directly relate 
to matter and satisfies the Laplace equation (e.g. see [8]). Also note that there is no Newtonian analogue to the 
magnetic Weyl tensor, which reflects the absence of gravitational waves within the limits of Newton's theory. 
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description, we need to supplement this set by an equal number of constraints. These come after 
contracting identity ([6b) with the permutation tensor of the space. Employing decomposition 
([2]), the result reads 

Scdad c a b d + d b u a - (d c co c )5 ab - i e abc 8 c e = . (12) 

The trace of the above, combined with the total antisymmetry of e abc , immediately leads to the 
familiar vorticity constraint 

d a uj a = , (13) 

guaranteeing that oo a is a solenoidal vector. On the other hand, taking the symmetric and 
trace-free component of Eq. (|12p we arrive at 

curkr afe + d( b u a) = , (14) 

where curlT a 6 = £ c d{ad c T b ) d for every symmetric and trace-free tensor of rank two. Finally, the 
antisymmetric part of (I12D leads to 

-d a @-d b a ab + cmhj a = 0. (15) 

and provides a relation between the gradients of the three kinematic variables. The reader is 
referred to Eqs. (1.3.6)-(1.3.8) in [3J for a comparison between the Newtonian and the relativistic 
kinematic constraints. Here, we simply note that in relativity oj a is not generally a solenoidal 
vector. Further discussion on Newtonian covariant hydrodynamics can be found in [3J. 

So far our analysis applies to all situations where the fluid description is valid. Typical 
cosmological models, for example, have > and u, a, p ~ 0. A non-rotating star, on the 
other hand, is characterised by w ^ and by 0, a ~ 0, while p oc p 7 (with 7 = constant) is 
a commonly used equation of state for the matter. Variations of the latter are also used in 
galactic studies, where observations indicate ~ and we can use Oort's constants to estimate 
the associated shear and vorticity. 

2.3 Perfect fluids 

When the gravitational field is specified and an equation of state for the fluid has been intro- 
duced, expressions ([5]), (|8l)- (fTT1) and (fT3]) - ([T5]) provide the fully nonlinear covariant equations 
that monitor the hydrodynamic behaviour of a self-gravitating Newtonian fluid. Here, we will 
consider the case of a barotropic perfect fluid with p = p(p) and 7r ab = 0. Under these conditions 
the Navier-Stokes equation reduces to 

A a = -^A a , (16) 
a 

where c^ = dp/dp is the square of the adiabatic sound speed and A a = (a/p)d a p. The latter is a 
dimensionless quantity that describes spatial variations (inhomogeneities) in the density of the 
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fluid, as measured between two neighbouring flow lines (e.g. see [6]). Assuming, for simplicity, 
that both the sound speed and the scale factor have zero spatial dependence, the above leads to 

d b A a = -^A ab , (17) 
cr 

with A ab = ad b A a . This variable is also dimensionless and, in contrast with the relativistic 
case, has zero skew part (i.e. Auy = 0). Thus, within the Newtonian framework, A ab can 
be used to describe density perturbations (by means of the scalar A = A a a = ad a A a ) and 
shape distortions (via the symmetric and trace- free tensor A^ = ad/j,A a \) but no vortex-like 
(i.e. vector) inhomogeneities. On using result {JTTJ) , expressions (f8))- (fT0j) take the form 

Q = - l -Q 2 -\np-^A-2{a 2 -u J 2 ) + K, (18) 
3 2 a z 

2 c^ 

&ab = — r ©Cafe - E ab 1 A< o6 ) - cr c ( a cr c fe ) + w c ( a a; c 6 ) (19) 



and 
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Wat = ~g @Uab ~ 2a c i a uj c b ] , (20) 

respectively. According to (|18p . overdensities (i.e. perturbations with A > 0) tend to enhance the 
gravitational collapse of the fluid, while underdensities support against it. In addition, following 
(fTUj) and ([2Pj) , the barotropic fluid can act as a source of shear anisotropy but does not generate 
vorticity. The same behaviour has also been seen in the relativistic studies (e.g. see [3]). Here, 
the main difference is that rotation remains unaffected by the fluid pressure (compare expression 
(|20p to Eq. (3.2.8) in [3]). As a result of this, which is due to the zero curvature of the Euclidean 
space, vorticity can never grow in expanding Newtonian models with vanishing shearH 

Finally, we note that one may monitor the acceleration or deceleration of an expanding 
Newtonian (barotropic) fluid by recasting Eq. (I18p into the form 

I 6\ = \k P + 4 A + 2(a 2 -co 2 ) -A, (21) 
3 2 a z 

where q = —ad/a 2 is the deceleration parameter. The above also shows how "voids", namely 
underdense regions with A < 0, tend to accelerate the expansion by acting together with the 
vorticity and the (positive) cosmological constant. 



2.4 Imperfect fluids 

One may look at the implications of fluid viscosity by considering an imperfect medium with 
nonzero anisotropic pressure. Maintaining the p = p(p) assumption of the previous section for 
simplicity, we introduce the phenomenological expression 

T^ab = -Acr a fe , (22) 

6 In general relativity, the rotational behaviour of the fluid also depends on its pressure. In particular, vorticity 
grows when the (dimensionless) adiabatic sound speed is greater than -y/2/3 (see [9] and also [3]). 
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with A = \(p,p) > being the viscosity coefficient (e.g see [1]). When the latter is a slowly 
varying function, the above combines with constraint (|15p to recast the momentum conservation 
law (see Eq. (03 )) into 

c? . A / 2 



4. = -— A B +- — 2a + curL; a , (23) 
a p \6a ) 

where Z a = ad a @ describes inhomogeneities in the volume expansion/contraction. Thus, by 
exploiting the advantages of the covariant expressions (in particular by involving constraint 
(|15p ). we were able to recast the viscosity term of ([5b) into a kinematical one. Proceeding as 
with the perfect fluid, we assume that both the sound speed and the scale factor depend solely 
on time. This allows the direct comparison of the two cases and leads to 

d b A a = A a6 + (Z ab - Z a A b ) + - (d 6 curla; a - - A fc curkj a ) , (24) 
a z 3arp p \ a J 

with Z ab = ad b Z a . Substituting the trace, the symmetric trace- free part and the skew component 
of the above into Eqs. (f8|)- (fT0|) . we arrive at 



1 _ 2 1 c? . 2A 



e = --e z - -Kp- -4 a + 



Z - ( Z a + y curl(j a ) A" 



3 2 y a? 3a 2 p 
-2(cr 2 - J 2 ) + A, (25) 



and 



Wa6 



2 



A_ 

ap 



2A 



# A<a6> + 3^V ^ <a6> ~ Z(aAb) 



(A( a curlw b) - a5( a curla; b )) - a c{a a c b) + uj c ^ a u c b) 



2 2A 
"3 6 ^-3^ 



^[ a A fe ] - - (aA [a curlw fc ] - a 2 a [a curla; 6 ]) 



(26) 



(27) 



respectively. Not surprisingly, we find that viscosity can modify every apsect of the model's 
kinematics in a variety of ways. Perhaps the most direct effect, relative to the barotropic- fluid 
case, is seen in Eq. (I27p . The latter shows that viscosity, together with an overall inhomogeneity, 
can act as a source of rotation (at the second perturbative level). 



2.5 Hydrodynamic flows as purely gravitational motions 

Keplerian motions are central to mass measurements. The observational determination of the 
masses of various astrophysical systems is usually based on the assumption of purely gravita- 
tional motions. For example, the central mass concentration in various galaxies is estimated by 
Doppler-shift measurements of radiative sources, which are assumed to move along Keplerian 
trajectories (e.g. see |10|). Nevertheless, there are known cases where the non-gravitational forces 
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are strong enough to affect these trajectories and where a hydrodynamic description of the mo- 
tion is more appropriate [11] . Then, one would like to know whether the standard measurements 
have overestimated or underestimated the available amount of matter. 

One way of addressing this question is by rewriting key hydrodynamic equations into a 
"Keplerian" form and then examining the implications of such a transformation for the dynamics 
of the physical system under consideration. Following [12J, and in absence of anisotropic pressure, 
we may combine Eqs. ([3]) and ([5b) to 

v a = -d a $ --dap. (28) 

p 

Setting V = 1/p as the specific volume, we introduce an equation of state of the form E = 
E(p, V), where E is the specific internal energy of the fluid. We may also define the associated 
temperature T = T(p, V) and specific entropy S = S(p, V) by (e.g. see [I]) 

dE + pdV = TdS. (29) 

where the right-hand side vanishes when adiabaticity holds. In the case of purely isentropic 
motions (i.e. when S is spatially and temporally constant), one can use the above expression to 
recast ([28]) as 

v a = -da$-d a (E+?j , (30) 

thus incorporating the internal properties of the fluid into Euler's equation. This means that 
isentropic hydrodynamic flows can be seen as entirely gravitational motions under the new, 
effective potential 

$ = <S> + E + -, (31) 
P 

which is shown to correspond to an effective mass-density given by 

d 2 $ = ± K p. (32) 

The "Keplerial" density introduced above can be expressed in terms of the wider fluid charac- 
teristics, like its internal energy and pressure, through definition (|3ip . In general, p is different 
from its hydrodynamic counterpart and their difference 

l - K {p-p) = d 2 {E + ^j , (33) 

depends on the aforementioned physical properties of the fluidQ This result also offers a way 
of measuring the "error-bars" between mass estimates based on purely gravitational motions, 
relative to those using the more realistic hydrodynamic approximation. For instance, if the 
"virtual" density p is smaller than the "actual" one (p), mass measurements using Keplerian 
motions will underestimate the available amount of matter. Although the results generally 
depend on the particulars of the physical system under consideration, there seem to exist realistic 
astrophysical environments where p < p (see [12[ [13] for further astrophysical discussion) . 

7 The effective mass density p does not generally obey a continuity equation of the simple form 
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3 Covariant magnetohydrodynamics 



Covariant techniques were introduced to the study of electromagnetic fields in [TJ] and more 
recently in [Taj (see also [TB] for an up to date review). All these studies are relativistic, however, 
and so far the Newtonian version of 1+3 covariant electrodynamics and magnetohydrodynamics 
(MHD) has been missing from the literature. 



3.1 Maxwell's equations 

In a two-fluid plasma description the charge carriers are the positive ions and the electrons, 
which are treated as two coupled conducting fluids. The matter density, the charge density and 
the current density of the one-fluid description are 

p = m + n + + m_n_ , q = e(n + — n_) (34) 

and 

J a = e(n + v+ - n-V~) , (35) 

respectively (e.g. see |17j). In the above e is the electron charge, m±, are the ion and the electron 
masses, n± represent their number densities and are the associated velocities. In the case of 
global electric neutrality we have n + = ri- and the centre of mass of the ion-electron system 
has the "bulk" velocity | 

v a = ■ {m + vt + m_u~) . (36) 

m + + m_ 

Within the single fluid approach and at the limit of resistive magnetohydrodynamics (MHD), 
the displacement current (dtE a ) is negligible. Then, Maxwell's equations reduce into a set of 
one propagation equation 

d t B a = -cnvlEa , (37) 

and three constraints 



curl5 a = J a , d a E a = and d a B a = , (38) 

having adopted the Heaviside-Lorentz electromagnetic units. The above, which respectively 
correspond to Faraday's law, Ampere's law, Coulomb's law and Gauss' law, are supplemented 
by Ohm's law. For a fluid with nonzero electrical resistivity, the latter reads 

J a = q(E a + e abc v b B c ) , (39) 

with <^ representing the (scalar) electrical conductivity of the medium. This form of Ohm's 
law corresponds to the resistive MHD approximation, which applies to fluids with small but 
finite electrical resistivity. In general, Eq. (139j) contains several additional terms - like those 
representing the Hall and the Biermann-battery effects (e.g. see expression (3.5.9) in [19J)@ 

8 See [18] for a generalisation to general relativity and further discussion. 

9 The potential implications of a non-conventional (anomalous) form of electrical resistivity were discussed 
in [20]. Also, for a comparison with the fully relativistic counterpart of (|39[) . the reader is referred to [21] . 
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Solving ([39]) for the electric field vector, substituting the result into Eqs. ([37]) . using decom- 
position ([2]), constraint (|38a ) and involving the convective derivative operator, we obtain the 
covariant form of the Newtonian magnetic induction equation 

B a = ~QB a + (a ab + uj ab )B b + -d 2 B a . (40) 

at the resistive MHD limit. Comparing the above to its relativistic counterpart (see Eq. (3.2.4) 
in [16]) we notice that the relative motion terms (i.e. the first two terms in the right-hand side of 
the two formulae) are identical. We also note the absence of the acceleration terms from (138a ) 
and (|40p - compare the former to (3.2.3) in [16]. This absence reflects the fact that Newtonian 
physics treats time and space as entirely separate entities. 

Similarly, employing ([39]) together with Ampere's law (see Eq. (f38a )). expression ([38b ) - 
Coulomb's law - assumes the covariant form 

2B a ta a = -v a J a , (41) 

suggesting that the sum v a J a = v a cuv\B a acts as an effective charge density relative to a rotating 
observer (when B a to a ^ - see also Eq. (3.2.5) in |16|). 

3.2 Magnetic evolution 

Contracting the magnetic induction equation (see (|40p ) along the field vector leads to the non- 
linear evolution formula of the magnetic pressure, namely to 



(b 2 \ = --&B 2 -2a ab U ab + -d 2 B 2 2 
V / 3 c 

where 



3 ? ? 



d{b B a)) -[d[b B a] 



(42) 



n a b = —B( a B b } = - B 2 5 ab — B a B b , (43) 

by definition. The latter is a symmetric, trace-free tensor that describes the magnetic anisotropic 
stresses and corresponds precisely to its relativistic counterpart. By definition, IL ab = M. ab — 
(B 2 /6)S ab , where M ab = (B 2 /2)5 ab - B a B b is the Maxwell tensor [23]. Thus, in agreement with 
the relativistic analysis (see section § 5.1 in [H]), the -B-field exerts an isotropic pressure equal 
to pb = M a a /3 = B 2 /6 and has an anisotropic pressure component given by U ab - We also 
note the quantities d^B^ and d[ b B a ] in the right-hand side of (j4"2|) . These may be respectively 
interpreted a the shear and the vorticity analogues of the -B-field - see Eq. (I55p below - and are 
important in highly distorted and turbulent magnetic configurations. 

Definition (|43p immediately ensures that U ab B b = -(2B 2 /3)B a . This in turn means that the 
-B-field is an eigenvector of the II a fe-tensor, with —2B 2 /3 being the associated eigenvalue. The 
negative sign shows that the magnetic pressure in the direction of the field lines is negative and 
reflects the tension properties of the latter (see also § !4.1l below). Projecting (|43p orthogonal to 
the magnetic forcelines, on the other hand, we find a positive eigenvalue equal to B 2 /3, which 
verifies that the field exerts a positive pressure in that plane [22 1 . In other words, every single 
field line acts like an elastic rubber band under tension, while neighbouring lines tend to push 
each other apart [23J. 
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Finally, following (|43p . it becomes immediately clear that the magnetic induction equation 
- together with expression (|42p - also monitors the time evolution of the anisotropic pressure of 
the field. On the other hand, the divergence of (|43|) provides the associated constraint, namely 

d b U ab = e abc B b cm\B c - \ 3 a B 2 . (44) 

6 

4 Resistive magnetohydrodynamics 

The resistive (or real) MHD scheme is believed to provide a good approximation to a variety of 
"typical" physical environments. For example, when the Larmor frequency and bulk velocity of 
the plasma are small, or when the dimensions of the system under study are large. 

4.1 The Lorentz force 

In covariant terms, the evolution of a non-relativistic magnetised plasma of finite electrical 
resistivity is monitored by the nonlinear set 

P=-&P, (45) 

P A a = -8 aP - 8\ ab - e abc B b cm\B c , (46) 

d 2 $ = ±K P , (47) 

consisting of the continuity equation, the Navier-Stokes equation and Poisson's formula respec- 
tively. These are supplemented by Maxwell's equations, which applied to an electrically resistive 
medium and written in covariant form read 

B a = -^QB a + (a ab + u; ab )B b + -d 2 B a , (48) 

curlB a = J a , (49) 

d a B a = , (50) 

and 

2B a tu a = -v a J a = -v a cuvW a . (51) 



The momentum conservation is reflected in ([46 p . which is the Navier-Stokes equation generalised 
to a (globally neutral) magnetised fluid. This expression can be obtained directly from its 
hydrodynamic counterpart (see ((3b)) by implementing the aforementioned fluid description of 
the £?-field. To be precise, Eq. (|46p emerges after replacing p with p + B 2 /Q and ir ab with 
■K ab + U ab in the right-hand side of |3b), while using constraint ((Hj) at the same time. Note 
that there is no magnetic contribution to the total inertial mass in the left-hand side of (|46p . 
or to the total gravitational mass in the right-hand side of Eq. (|55|) - see § 14.21 below. This is 
a significant change, with respect to the relativistic case (compare to expressions (5.3.3) and 
(5.5.1) of |16j), which implies that there is no Newtonian analogue to the relativistic energy 
density of the 5-field. Finally, following (|46p . we note that the sum A a B a has no magnetic 
dependence. This ensures that the magnetic field has no effect along its own direction. 
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The set (j4"5j) - ([5Tj) is supplemented by the kinematic propagation and constraint equations 
(f8j)- (fl"5j) . once the latter have been appropriately adapted to our electrically resistive magnetised 
environment. Within the limits of the Newtonian theory, the above named formulae contain no 
explicit magnetic terms. This means that the kinematic effects of the -B-field propagate solely 
through the fluid acceleration and specifically via the Lorentz-force term in the right-hand side 
of Eq. (|46p f^l For a globally neutral medium, the Lorentz force depends exclusively on the 
-B-field and splits into two stresses according to 

e abc B b cmW c = ± 8 a B 2 - B%B a , (52) 

where the first term in the right-hand side is due to the isotropic pressure of the field (see 
§ 13.21 below) and the second carries the effects of the magnetic tension. The tension stress also 
reflects the elasticity of the field lines and their tendency to remain straight [23J. When these 
two stresses balance each other out, the -B-field reaches equilibrium. 



4.2 Nolinear kinematics 

Proceeding as in § 12.31 we ignore the anisotropic pressure of the fluid and assume a barotropic 
medium by setting p = p(p). Then, the MHD version of the Navier-Stokes equation (see (146]) ) 
takes the form 

A a = -^A a -^B a + -B b 8 b B a , (53) 
a 2a p 

with c 2 = B 2 /p and B a = (a/ B 2 )d a B 2 . The former is the Alfven speed, which determines the 
propagation of MHD disturbances and also provides a measure of the relative strength of the 
-B-field. The latter is a dimensionless variable that describes spatial variations in the (isotropic) 
magnetic pressure. Note the last term in the right-hand side of Eq. (|53p . which carries the 
effects of the magnetic tension (see decomposition (J52j) above). When the sound speed, the 
Alfven speed and the scale-factor have a spatially homogeneous distribution, the gradient of the 
above leads to 

d b A a = - 4 A ab - ^% B ab - — A b B c 8 c B a + - d b B c 8 c B a + - B c d c 8 b B a , (54) 
a z zar ap p p 

where B ab = ad b B a . The overall magnetic effect is rather involved and propagates via the last 
four terms. Of these, the first is triggered by the isotropic pressure of the field and the rest are 
due to the tension properties of the magnetic forcelines. 

Substituting the trace of (|54p into Eq. ([8]), we obtain the nonlinear form of Raychauduri's 
formula for a magnetised, self-gravitating Newtonian fluid of zero total charge. In particular, 
using constraint (138b ). we arrive at 

9 = -\® 2 -\ K p- C \A-^B--A*B b d b B a -2(o 2 -a 2 B ) 
3 2 or 2a £ ap \ ' 

+2 (u? - lo 2 b ) , (55) 

10 The lack of explicit magnetic terms in the propagation formulae <JHJ _ (tH]) an d the absence of acceleration terms 
in Eqs. (|13p - (|15p . represents a considerable change relative to the relativistic case (see [16, for details). 
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where B = B a a , o\ = d {b B a) d^ b B a ^ /2p and w| = d [b B a] d^ b B a y2p. The former describes scalar 
variations in the magnetic pressure, while the last two can be interpreted as the magnetic 
analogues of the shear and the vorticity respectively. According to the above, the compression 
of the field lines (which corresponds to an increase in the magnetic pressure and B > 0) assists the 
gravitational pull of the matter. The dilution of the magnetic forcelines, on the other hand, acts 
against contraction. We also note that the effect of the magnetic shear and vorticity opposes that 
of their kinematic counterparts (see also [21]). The reason behind this counterintuitive behaviour 
is the magnetic tension. Both erg an d ujb are triggered by the elasticity of the magnetic forcelines 
and therefore react to any agent that distorts them. The magnetic vorticity, in particular, is 
the response of the field's tension to the twisting of its forcelines. The resulting stress slows 
the rotation down and this effect is commonly referred to as "magnetic braking". Analogous 
behaviour has also been observed in relativistic studies (see [25J for more details and further 
discussion). The key difference here, as a result of the Euclidean nature of the Newtonian space, 
is the absence of the general relativistic magneto-curvature stresses. 

Substituting the symmetric and trace-free component of the auxiliary expression (|54p into 
the right-hand side of ([9]), leads to 

2 c^ c 2 1 1 

Pah = —^&^ab-E ab -^A^--^B^ + -B c d c d^B a ^ + -d^ b B c d c B a ' } 

-— B c A (a d c B b) - a c{a a c b) + uj c{a uj c b) . (56) 

The above shows how anisotropies in the distribution of the magnetic pressure and in that of 
the field gradients affect the evolution of the kinematic shear. In particular, despite the lack of 
a direct contribution from the magnetic anisotropic pressure, the -B-field acts as a shear source 
in a variety of ways 1^1 Note that of the four magnetic source-terms in Eq. (|56p . the first is due 
to the field's pressure and the last three are the result of its tension. 

Finally, the skew part of decomposition (|54f) . together with the (strictly Newtonian) results 
B[ab] = = Ar ob i, transform Eq. (fTT|) into 

wa = 4 0w «-r B b d b curlB a - ±- e abc d b B d d d B c + e abc B d A b d d B c + a ab co b . (57) 
6 2p 2p 2ap 

This expression reveals the role of the -B-field as a source of rotation, either on its own or through 
its coupling to the density gradients. It should also be noted that there are no effects due to the 
isotropic magnetic pressure in Eq. (|57p . with all the i?-terms coming from the field's tension. 
Following (I56p and (157j) . even if the fluid is originally shear-free and non-rotating, it will not 
remain so once a magnetic field is introduced. 

We close this section by noting that, according to Eqs. (I13j) -(|15j). the kinematic constraints 
contain no explicit magnetic terms and therefore are only indirectly affected by the field's pres- 
ence. We should also underline the benefits from using the covariant approach. These are 
multiple because the formalism streamlines the equations, while maintaining maximum detail 
and physical transparency. Finally, we note that the expressions given in § H] can be used to 



1 Recall that in relativistic studies the magnetic iI a i,-tensor is an explicit source of shear anisotropies [16] 
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study the Newtonian evolution of any electrically resistive and globally neutral fluid in the pres- 
ence of a magnetic field. In addition, the formalism developed so far can be extended to study 
the behaviour of inhomogeneities, both at the linear and at the nonlinear level. 

4.3 Hydrodynamic reduction of magnetised flows 

Ideal fluids have zero anisotropic pressure by definition. When, in addition, the tension compo- 
nent of the magnetic Lorentz force is also zero (i.e. for B b d b B a = - see decomposition (|52|) ). 
the generalised Navier-Stokes equation simplifies to 

pA a = -8 aP - 1 8 a B 2 . (58) 

Realistically speaking, the above is only an approximation and holds when the Lorentz force 
is dominated by the (positive) pressure of the 5-field. In such an environment, the non- 
gravitational acceleration of the fluid (i.e. the vector A a ) comes purely from a potential. Then, 
the MHD motion reduces to a simple hydrodynamic flow with 

pA a = -8 a P , (59) 

where the scalar P = p + B 2 /2 acts as an effective hydrodynamic pressure (see Eq. (I53p ). This 
new motion is monitored by the formulae of § 12.21 after replacing expression ([5b ) with (|59p 
and the pressure of the original fluid with the above given effective pressure P. One can also 
go a step further and use the transformations of § 12.51 to represent the MHD flow of (|58p , (|59p 
as a "purely gravitational" motion. This time the effective potential will also depend on the 
magnetic pressure. 



5 Ideal magnetohydrodynamics 

In most astrophysical and cosmological studies, magnetic fields are treated within the limits of 
the ideal MHD approximation. The latter, applies to highly conductive media with essentially 
zero electrical resistivity. Although overly idealised and simplistic, the perfect MHD scheme still 
seems to provide the correct description in a variety of studies. 

5.1 Maxwell's equations 

When dealing with a perfectly conductive medium, namely at the <j — > oo limit, the Ohmic 
current in Eq. (|39|) vanishes (i.e. J a /s — > 0) and the associated electric field is given by the 
simple expression 

E a = -e abc v b B c . (60) 

In these environments the nonlinear equations monitoring a globally neutral, self-gravitating 
Newtonian fluid in the presence of a magnetic field are 

P = -&P, (61) 
pA a = -d aP - d\ ab - e abc B b cm\B c , (62) 
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d 2 $ = -K P , (63) 

B a = - ~ QB a + (a ab + co ab ) B b , (64) 

curlB a = J a , (65) 

d a B a = , (66) 

2B a u a = -v a J a = -v a cui\B a . (67) 

Relative to the resistive-MHD case of the section § HI we note the absence of a diffusion term 
in the right-hand side of the induction equation (compare Eqs. (|40p . (|48p to expression (|64p 
above) l 12 l This guarantees that the magnetic field lines remain frozen-in with the matter. In 
particular, (|64p ensures that X a = a?B a is a relative position vector connecting the same particles 
at all times (i.e. X a = X b d b v a - see footnote [3] in § 12.11 and also [6| fTB]). 



5.2 Magnetic evolution 

Relation (|64p also shows that, in the absence of shear anisotropies, the magnetic strength either 
dilutes with the expansion or increases with the contraction of the fluid. Then, recalling that 
0/3 = a I a, the magnetic induction equation reduces to 

B a = -2 Q B a . (68) 

An immediate consequence ia that the magnetic flux, here represented by the quantity a 2 B a , 
remains conserved in time. Moreover, the ideal-MHD counterpart of Eq. (|42h reads 

(B 2 )' = -^QB 2 -2a ab U ab , (69) 

with U ab given in (|43p . Therefore, for zero shear anisotropy, we recover the familiar from 
the relativistic studies radiation-like evolution (i.e. B 2 oc a -4 ) of the magnetic pressure. The 
presence of shear, on the other hand, will generally modify the aforementioned "adiabatic" 
pattern. This can happen during the realistic (i.e. anisotropic) collapse of a magnetised proto- 
galactic cloud and lead to the amplification of the embedded -B-field beyond the limits of the 
simple spherical-collapse models [26J. 

As mentioned in § 13.21 the magnetic induction equation monitors the time evolution of both 
the isotropic and the anisotropic pressure of the 5-field. At the ideal-MHD limit the time 
derivative of f|43[) combines with expressions (|68p and (|69p to give 

1U = ~ en ab + 2n c{a a c b) - 2U c{a u; a b) - I B 2 a ab , (70) 

while the associated constraint is still given by (|44p . In the absence of shear and vorticity, the 
above leads to H ab oc a -4 , in line with the evolution of its isotropic counterpart. Thus, when 
the anisotropy is small, the i?-field has a radiation-like evolution to first approximation. 

12 The kinematics of a perfectly conductive ideal fluid are still monitored by the "resistive" formulae of § 14.21 
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Turning to the kinematics of perfectly conducting media, we note that the magnetic effects 
on a (globally neutral) fluid propagate via the Lorentz-force term in the right-hand side of 
the generalised Navier-Stokes formula (see Eq. (|62|) i. The form of the latter is independent 
of the electrical resistivity of the matter, since it contains no related terms. This means that 
relations (|55p -(|57|). together with constraints (|13|) - (|15| ). also govern the kinematics of an ideal- 
MHD medium] 13 ! When an equation of the state for the matter is introduced, these expressions 
monitor the nonlinear evolution of the magnetised medium completely and in a fully covariant 
manner. 

We finally note that, when the fluid is perfect, the magnetic field becomes the sole source 
of anisotropy. The magnetically induced effective viscosity can be related to that of the shear 
in a way that closely resembles the phenomenological equation of state introduced in § 12.41 (see 
Eq. (|22p there). Thus, assuming that the -B-field is a shear eigenvector, we may set (J a \ ) B b = 
(2fj,/3)B a , where 2/i/3 is the associated eigenvalue. Also, following definition (|43p . we find that 
Ii ab B b = -(2B 2 /3)B a and subsequently arrive at 

Ii ah = -Xa ab , (71) 

with A = B 2 /fi acting as an effective coefficient of magnetic viscosity |22j. 

6 Discussion 

Newtonian theory offers a very good approximation to general relativity in weak-gravity envi- 
ronments and also on scales well inside the Hubble radius. The covariant approach to Newtonian 
hydrodynamics is a Lagrangian description based on a relative-motion treatment that exploits 
the irreducible kinematical quantities of the motion. Although the formalism was originally 
applied within the framework of Newton's theory, it has since been used primarily in relativistic 
studies. On the other hand, while the 1+3-covariant techniques have been employed for the 
study of relativistic electromagnetic fields, so far a Newtonian version of that work has been 
missing. 

The present paper reviews and extends the existing work on Newtonian covariant hydrody- 
namics on the one hand, while on the other it applies the covariant techniques to magnetohy- 
drodynamic studies. Exploiting the advantages of the relative-motion treatment, we supplement 
the standard hydrodynamic formulae with a set of three propagation and three constraint equa- 
tions that monitor the evolution of the irreducible kinematical variables. The latter, namely 
the volume expansion/contr action, the shear and the vorticity, describe the relative motion of 
neighbouring flow lines. The aforementioned formulae are obtained in a manner analogous to 
that of their relativistic counterparts and this facilitates the direct comparison of the two sets. 
In fact, the close analogy between the Newtonian and the relativistic equations is maintained 
throughout the paper and this allows the unambiguous identification of their differences. 

We consider perfect as well as imperfect (viscous) media, looking for differences in their 
kinematical behaviour. Not surprisingly, the extra degree of freedom that viscosity introduces, 

13 Directly, the electrical resistivity of the medium affects only the evolution of the embedded .B-field (compare 
Eqs. (|48[) and ((64} ) . The latter then carries these effects to the kinematics and the dynamics of the magnetised 
medium. 
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means that the kinematics of a viscous fluid are considerably more involved. Following our 
analysis, the key contribution of viscosity is perhaps through its role as a source of rotation. 
Focusing on isentropic fluids, we also discuss how hydrodynamic flows can be represented as 
"purely gravitational" motions due to a new (effective) potential. The latter incorporates ad- 
ditional characteristics of the fluid, like its internal energy and pressure, and corresponds to a 
new (effective) mass density. The relation between the "actual" (the hydrodynamic) and the 
"virtual" (the effective) mass density has been used to estimate the accuracy of astrophysical 
mass measurements based on the assumption of purely gravitational (Keplerian) motions. 

Assuming an imperfect MHD fluid of zero total charge, we derive the covariant version of 
Maxwell's equations within the limits of Newtonian gravity. In an environment of small but 
finite electrical resistivity, we monitor the evolution of the magnetic field completely. This 
means providing the nonlinear propagation and constraint equations for both the isotropic and 
the anisotropic magnetic pressure. The compactness of the covariant formalism also allows us to 
identify the impact of the -B-field on the kinematics of the fluid in detail. In practice, this means 
isolating the effects due to the ordinary (the positive) magnetic pressure, from those coming 
from the tension of its forcelines. The former affect the volume evolution and also the shape of 
a given fluid element, but not its rotational behaviour. The impact of the magnetic tension, on 
the other hand, is more widespread and sometimes counterintuitive. Thus, the elastic properties 
of the field lines are shown to act as sources of rotation, either on their own or through their 
coupling to density inhomogeneities. In an analogous way, magnetism is also found to trigger 
shear distortions. Moreover, when looking into the implications of the -B-field for the volume 
evolution of the fluid, we identify magnetic analogues of the shear and the vorticity. Both carry 
the tension properties of the field and oppose the effects of their kinematic counterparts. The 
"magnetic vorticity" term, in particular, tends to slow the rotation down and leads to what 
is commonly referred to as magnetic braking. Finally, we apply our analysis to the limit of 
ideal magnetohydrodynamics and also discuss how certain MHD flows can be reduced to simple 
hydrodynamic ones. 

The formalism developed here can be applied to a variety of astrophysical and cosmolog- 
ical environments, where the Newtonian theory is a good approximation. This includes non- 
relativistic astrophysical MHD and galaxy formation studies. In the latter case, for example, 
one could use the linearised version of our equations to follow the linear regime of a magnetised 
protogalactic cloud (with size well below the horizon scale) . Similarly, the full expressions can be 
employed to monitor the nonlinear evolution of the protogalaxy, once the latter has decoupled 
from the background expansion and started collapsing. More specifically, our equations will 
enable one to look for effects outside the limits of the ideal MHD. The latter are expected to 
play a role during the nonlinear regime of galaxy formation (at least locally). Applications of 
this type will be the subject of future work. 
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